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The question of determining the maximal number of mutually unbiased bases in dimension six 
has received much attention since their introduction to quantum information theory, but a definitive 
answer has still not been found. In this paper we move away from the traditional analytic approach 
and use a numerical approach to attempt to determine this number. We numerically minimise a 
non-negative function Did, at of a set of JV+1 orthonormal bases in dimension d which only evaluates 
to zero if the bases are mutually unbiased. As a result we find strong evidence that (as has been 
conjectured elsewhere) there are no more than three mutually unbiased bases in dimension six. 



I. INTRODUCTION 



Mutually unbiased bases have uses in a variety of top- 
ics in quantum mechanics. The first person to consider 
their use was Ivanovic [IJ who introduced them in the 
problem of state determination. Here it is found that 
measurements based on mutually unbiased bases are opti- 
mal for constructing the state of an ensemble of systems. 
This application was investigated further by Wootters 
and Fields some time later 0). In 2002 Cerf, Bourennane, 
Karlsson and Gisin Q explicitly use mutually unbiased 
bases to extend the BB84 [H and the six-state proto- 
col Q quantum cryptographic schemes (25[. The Mean 
King 's Problem is also related to mutually unbiased bases 



Mutually unbiased bases are simple to define. Sup- 
pose we have a c?-dimensional Hilbert space and a number 
of orthonormal bases for this space described by vectors 
\ipkm): where m labels one of the vectors in basis k. We 
call these bases mutually unbiased if 
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\(lpk,m\lpl,n)\ = ~/= 



Vfc 7^ I, Vm, n = 1, . . . , d. (1) 



A commonly asked question is then: What is the max- 
imum number of mutually unbiased bases that can be 
found for general dl 

Ivanovic, in his paper on state determination con- 
structed a set of d + 1 mutually unbiased bases for all 
prime d. Wootters and Fields went further and proved 
not only that d+1 is the maximal possible number of mu- 
tually unbiased bases in dimension d, but that this bound 
is realised for all prime power dimensions. 

The smallest dimension for which the maximum num- 
ber of mutually unbiased bases is unknown is d = 6. It is 
known Q that for any dimension a the maximal number 
of mutually unbiased bases 97t(d) can be bounded below 



e.g. 



97t(d) >min(97t(^ ),..., Tl( P e /)) 
= min(p? + l,...,p? + l), 



where d = p^ 1 . . . p e r T is the prime factorisation of d. Thus 
it is clear that 971(6) > 3. A number of the constructions 
for three mutually unbiased bases in d = 6 have been ex- 
amined by Grassl and it has be en p roven that these 
cannot be extended further. Archer [ljj also proved that 
the methods used to determine 971(d) for prime power 
dimensions cannot be generalised to non-prime power di- 
mensions. Bengtsson et al. [l2j consider a constructions 
from a number of known types of Hadamard matrices but 
are unable to find more than three mutually unbiased 
bases for d = 6 from these. Other results establishing an 
improved lower bound for square dimensions exist [1 31 ] . 
but the general question is far from answered, and the 
above inequality is the only known bound for d = 6. 

The inability to analytically find four or more mutu- 
ally unbiased bases in dimension six has brought about 
a general belief that 971(6) = 3 (as conjectured in [HI). 
Rather than approach this problem analytically, we will 
reformulate it in such a way that we can perform a nu- 
merical minimisation to provide evidence to support or 
disprove this conjecture. This approach has been sug- 
gested by Bengtsson [l2j, and a similar calculation has 
been carried out independently in [l5[; however, to our 
knowledge, no literature presenting detailed results of a 
numerical investigation exists. 

The rest of this article is organised as follows. In sec- 
tion 2 we discuss our reformulation of the problem into a 
minimisation problem and the algorithm we use to find 
the minima. In section 3 we present the results of our 
various searches and analyse them. Finally in section 4 
we discuss other potential approaches to the problem and 
summarise our results. We also provide more technical 
details of the minimisation algorithm in the appendix. 



II. NUMERICAL METHODS FOR MUTUALLY 
UNBIASED BASES 
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In this section we show how the problem of finding mu- 
tually unbiased bases can be cast in the form of an opti- 
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misation problem. We will illustrate how this problem is 
constructed and how it can be attacked numerically. 



A. Constructing the optimisation problem 

Let us first state the general idea we are trying to pur- 
sue here. We wish to create a function of N + 1 orthonor- 
mal bases (if they existed) that would achieve a global 
minimum when these N + 1 bases are mutually unbiased. 
We then hope to obtain the minimum using an appropri- 
ate optimisation algorithm. 

We start with some useful preliminaries. Let \4> m ) = 



-\ ymn | 



be an orthonormal basis for 



,,d. 



Define a d x d matrix U by the coefficients of the indi- 
vidual vectors i.e. U nm = <p mn (so that column m corre- 
sponds to \(f> m ))', then the fact that the basis is orthonor- 
mal implies that U is unitary. Now suppose that A, B 
are two unitary matrices representing two orthonormal 

bases \a m ) = £) n =i A nm \n) and \(3 m ) = £) n =i B nm \ri) 
respectively. Then 



( A ^ B )mn = Y, A vmB P n = (a m \P n ) 



(2) 



then 



and so if the two bases are mutually unbiased, 
| (A^B) mn | = 1/Vd for all m, n = 1, . . . , d. 

So, we can represent N + 1 mutually unbiased bases 
by unitary matrices U±, . . . , Un, Un+i = 1 (where with- 
out loss of generality we have rotated every basis so that 
the last basis is simply the standard basis {\k)}f =1 ). So 
now we are looking for a function of the matrix elements 

( UfcUi ) that is minimised when each of these norms 
is equal to 1/Vd. One such function is 



m d , N (u u ...,u N ) = 

\<k<l<N+lm,n=l v 



(3) 



which has a minimum of zero when all of the matrices 
represent a set of N mutually unbiased bases [26[. So 
in principle we can determine whether a set of N + 1 
mutually unbiased bases exists by finding the global min- 
imum of this function over all sets of unitaries Uk- Using 
the fact that we can write any unitary matrix U in the 
form U = e iH , where H is Hermitian, we can state this 
minimisation in an unconstrained form i.e. 



can be solved using the Levcnbcrg-Marquadt algorithm 
[HI, EH- However, like many such algorithms, it is 
only guaranteed to find a local rather than a global max- 
imum, so to try and resolve this issue we will run the 
algorithm many times from a randomised starting point. 

The computational mathematics package MATLAB 
has an implementation of both this algorithm and a fast 
and accurate matrix exponentiation routine, and so we 
choose to use this package to implement our algorithm. 
Some of the more technical details of the implementation 
are given in the Appendix. 



III. RESULTS 

Our tests for d = 6 concentrated on both finding a full 
set of d+ 1 mutually unbiased bases, which exist in prime 
power dimensions, and also finding four MUBs which is a 
less computationally complex problem. For comparison 
we also conducted a number of tests for d = 2, 3, 4, 5, 7. 
In all tests we consider a set of mutually unbiased bases to 
be found if a set of unitary matrices are found such that 
%ld,N(Ui, . . . , Un) < 10~ 6 . The results are presented in 
Figure Q] 

Our attempts to find four mutually unbiased bases for 
d = 6 provide us with strong evidence towards their non- 
existence. Although the success rate was not high, our 
minimisation method is able to find four mutually un- 
biased bases for d = 4, 5, 7. The method was unable to 
find four mutually unbiased bases in dimension six and 
Figure [2] presents the frequency of minima obtained. The 
peaks suggest that in general our algorithm is converg- 
ing to a small number of local minima. Moreover, more 
than two-thirds of the test runs converged to the value 
0.051249, the minimum value of 0^6,3 our algorithm was 
able to find. 

For d < 6 the success rate of finding d + 1 mutually 
unbiased bases was very high. When attempting to find 
seven mutually unbiased bases for dimension six we ob- 
tain a similar picture as when trying to obtain four. In 
this case just under a third of the runs converged to our 
minimum value 1.584472. The long run time for minimis- 
ing 0^7^ [27l | prevents us from obtaining a larger number 
of results for d = 7, but on a small number of occasions, 
a full set of mutually unbiased bases is found. On in- 
specting the distribution of the obtained results, we do 
not observe the same grouping of the non-zero minima 
as in Figure [21 suggesting that the algorithm may not be 
converging to a local minimum in most cases. 



Minimise W diN {e lHl e l " N ) 



(4) 



B. Implementing the minimisation numerically 

The function 9T(i7i, . . . , Un) has the form of a sum of 
squares of a series of functions, and so the minimisation 
problem stated in Q is a non-linear least squares op- 
timisation problem. This kind of optimisation problem 



IV. CONCLUSION 

The numerical results we have presented here repre- 
sent clear evidence that a set of more than three mu- 
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d + 1 MUBs 


4 MUBs 


Dimension d 


No. of tests 


d + 1 MUBs found 


% success 


No. of tests 


4 MUBs found 


% success 


2 


2500 


2500 


100% 


N/A 


N/A 


N/A 


3 


2500 


2499 


99.9% 


N/A 


N/A 


N/A 


4 


2500 


2500 


100% 


2500 


2500 


100% 


5 


2500 


2495 


99.8% 


2500 


1510 


60.4% 


6 


2500 





0% 


10000 





% 


7 


250 


3 


1.2% 


3000 


26 


0.9 % 



FIG. 1: Results illustrating number of times each minimisation problem converges to zero, i.e. a set of mutually unbiased bases 
is obtained. 



tually unbiased bases do not exist in dimension six. It 
is also clear from the results in Figure [TJ that this ap- 
proach does not scale well at all as the dimension of the 
bases increases. However, we have presented only one 
possible approach to this problem. Our choice of min- 
imising function is by no means unique. Another plausi- 
ble test function could be based on the Shannon entropy 
H(xi,...,Xd) = -Y^^ilogcci, which, when £\ x t = 1, 
is maximised when Xi = 1 / d. We would of course have to 
use a different optimisation routine in this case. Another 
possibility is to simply paramcterise the unitary matri- 
ces Ui elementwise, and impose the unitarity condition 
as a constraint on the elements of each of the matrices. 
We have considered these approaches briefly, and found 
for the calculations we considered that our method was 
fastest. We also chose to approximate first derivatives 
via finite differences rather than calculate them directly; 
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FIG. 2: Histogram of numerically determined minimal values 
of 916,3- The separation of the histogram bars indicates that 
we really are obtaining local minima from our calculations. 



again, this was faster, and appeared to have no effect on 
accuracy, but this may not be the case with other meth- 
ods. However, given that we are working with a func- 
tion of 0(d 2 N) variables, minimising O^iv (or indeed 
any other appropriate function) efficiently and effectively 
may well become very difficult using any algorithm as d 
and N increase [28j. An example of a different method 
(although presented in less detail) can be found in [l5| 
(see also |29j). 

We can also consider the more general question of how 
close a set of N+l orthonormal bases can be to being mu- 
tually unbiased. The minimisation function used here is 
simply one way of quantifying how close a particular set 
of orthonormal bases is to being mutually unbiased (al- 
beit one that has a nice geometric interpretation [3.[l9|). 
However, there are clearly other ways of doing this. One 
possibility is to relax the actual definition of mutually 
unbiased e.g. we could define two vectors to be approxi- 
mately mutually unbiased if — e+ < 1(^10)1 — 1/V" < e+ 
for some e± > 0. In principle, one could determine what 
e± have to be to allow the existence of N + 1 approx- 
imately mutually unbiased bases \ipk,m) by considering 
the minimum and maximum value of \(ipk,m\'4'i,n}\ —1/d 
over all sets of N + 1 orthonormal bases. However, since 
these functions are not differentiable at many points, 
there may well be difficulties in finding the local extrema 
of these minimum and maximum functions (our brief at- 
tempt at trying this numerically yielded problems even 
for 3 bases in dimension 2). Knowing how close a set 
of orthonormal bases can get to being mutually unbiased 
could be helpful for example in determining efficient state 
determination protocols, as the information gained from 
a set of measurements in these bases simply depends on 
the geometrical relationships between the bases Q. 
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APPENDIX A: TECHNICAL DETAILS OF 
IMPLEMENTATION 

Our implementation as stated uses the MATLAB im- 
plementation of the Lcvenbcrg-Marquadt non- linear least 
squares minimisation routine (lsqnonlin). For speed, 
the evaluation of the function OLj.jv {U\, . . . , Un) was 
coded in C, and compiled within MATLAB. We set 
the routine to terminate when the value of the function 
^d,jv(^ii • ■ ■ i Un) changes by a value less than 10 -8 be- 
tween iterations (this threshold was chosen to give a de- 
cent balance between accuracy and time efficiency). 

The Levenberg-Marquadt algorithm only allows op- 
timisation over the reals, and so to deal with this we 
must parameterise the Hcrmitian matrices Hk in (jl]) 
using real variables. The C function that evaluates 



(4 



which are 



9^(2, at (EA, • ■ ■ 5 Un) then reconstructs the unitaries Uk and 
evaluates the real valued norms 
then passed back to the optimisation routine. 

As mentioned in the main body of the text, we pick a 
random starting point for each run of the algorithm. To 
generate the initial points, we must generate a number of 
random unitaries Uk , and then take the matrix logarithm 
of these unitaries to determine the values of Hk- We gen- 
erate random unitaries in the Haar measure (the analog 
of the uniform distribution for compact matrix groups, 
which can be implemented using the QR factorization. 
The following MATLAB code generates a random uni- 
tary matrix U in the Haar measure [20| : 

[Q,R] = qr (randn(dim) + i*randn(dim) ) ; 
U = Q*diag(exp(2*pi*i*randn(dim, 1) ) ) ; 

To exponentiate and take logarithms of matrices we 
utilise the in-built MATLAB functions. 
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